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State-of-charge  (SoC)  and  state-of-health  (SoH)  define  the  amount  of  charge  and  rated  capacity  loss  of  a 
battery,  respectively.  In  order  to  determine  these  two  measures,  open-circuit  voltage  (OCV)  and  internal 
resistance  of  the  battery  are  indispensable  parameters  that  are  obtained  with  difficulty  through  direct 
measurement.  The  motivation  of  this  study  is  to  develop  an  online,  simple,  training-free,  and  easily  imple- 
mentable  scheme  that  is  capable  of  estimating  such  parameters,  particularly  for  the  lithium-ion  battery 
in  battery-powered  vehicles.  Based  on  an  equivalent  circuit  model  (ECM),  the  electrical  performance  of  a 
battery  can  be  formulated  into  state-space  representation.  Also,  underdetermined  model  parameters  can 
be  arranged  to  appear  linearly  so  that  an  adaptive  control  approach  can  be  applied.  An  adaptation  algo¬ 
rithm  is  developed  by  exploiting  the  Lyapunov-stability  criteria.  The  OCV  and  internal  resistance  can  be 
extracted  exactly  without  limitations  of  a  system  input  signal,  such  as  persistent  excitation  (PE),  enhanc¬ 
ing  the  method  applicability  for  vehicular  power  systems.  In  this  study,  both  simulations  and  experiments 
are  established  to  verify  the  capability  and  effectiveness  of  the  proposed  estimation  scheme. 
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1.  Introduction 

Specific  energy  and  power  in  addition  to  cycle  life  are  pri¬ 
mary  concerns  of  batteries  utilized  in  electric  vehicles  (EV).  The 
lithium-ion  battery  is  believed  to  be  able  to  potentially  meet  these 
requirements  in  the  future.  Hence,  several  well-known  types  of 
lithium  batteries,  such  as  lithium  iron  phosphate,  lithium  poly¬ 
mer,  and  nano-phosphate  lithium-ion  cells,  have  been  developed. 
In  addition,  their  reliability  and  durability  are  very  susceptible  to 
operational  and  environmental  conditions,  particularly  in  EV  appli¬ 
cations.  To  achieve  the  required  EV  traction  power  and  range, 
low-voltage  lithium-ion  cells  are  generally  connected  in  series 
and  in  parallel  to  construct  a  dedicated  battery  pack.  A  battery 
management  system  (BMS),  along  with  protective  circuitry  and 
a  communication  bus,  is  provided  for  management,  monitoring, 
and  diagnosis.  Measurement  of  state-of-charge  (SoC)  is  one  of  basic 
functions  of  the  BMS,  which  indicates  the  remaining  charge  of  the 
battery  so  that  the  driver  can  be  reminded  to  charge  the  battery 
prior  to  its  depletion.  It  is  known  that  inaccurate  SoC  determina¬ 
tion  is  likely  to  cause  catastrophic  situations  such  as  overdischarge 
or  overcharge.  In  addition,  state-of-health  (SoH),  which  defines  the 
battery  performance  relative  to  its  fresh  condition,  is  used  to  predict 
its  end-of-life  and  aging.  In  another  application,  SoH  is  a  measure 
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used  to  analyze  the  effectiveness  of  the  battery  when  incorporated 
with  other  power  sources,  i.e.,  hybrid  energy  storage  systems,  in 
order  to  extend  battery  cycle  life.  Generally,  the  SoH  function  is 
uncommonly  seen  in  existing  BMS’s.  However,  there  is  an  increas¬ 
ing  need  to  monitor  the  battery  performance  in  EV  applications  so 
as  to  prevent  instant  breakdown  while  driving. 

Several  methods  for  SoC  determination  have  been  proposed 
since  1 938  [  1  ].  Such  methods  can  be  generally  categorized  as  direct 
or  indirect  approaches.  One  of  the  direct  approaches  simply  indi¬ 
cates  the  remaining  capacity  by  using  online  current  integration  [2  ]. 
Some  studies  have  improved  the  previous-demonstrated  accuracy 
by  compensating  for  temperature  and  aging  effects  [3,4].  On  the 
other  hand,  the  indirect  methods  determine  the  SoC  by  using  the 
battery’s  intrinsic  relationship  between  the  SoC  and  some  electri¬ 
cal  parameters  such  as  open-circuit  voltage  (OCV)  and  impedance. 
For  instance,  the  battery  OCV  naturally  declines  proportionately 
with  the  energy  expenditure  and  is  widely  used  for  SoC  indication 
[5-7].  Relying  on  an  impedance  spectroscopy  analysis,  the  SoC  and 
SoH  can  be  accurately  determined  [8,9].  In  recent  years,  the  con¬ 
cept  of  adaptive  estimation  based  on  neural  networks  [10,11]  or 
Kalman  filtering  [12,13]  has  been  proven  to  be  effective  in  SoC- 
and  SoH-determination  problems. 

The  SoH  indication  of  the  battery  has  attracted  growing  atten¬ 
tion  in  the  last  decade  for  automotive  applications,  e.g.,  EV  or 
hybrid  EV  (HEV).  Typically,  the  SoH  can  be  estimated  through  bat¬ 
tery  usage  history  or  performance  parameters.  The  usage  history 
might  involve  normal  charge-discharge  cycles  and  abusive  expe¬ 
riences,  e.g.,  instances  in  which  the  operational  limits  of  voltage, 
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current,  and  temperature  were  exceeded.  A  data  fusion  technique 
is  utilized  to  integrate  the  aforementioned  data  for  the  estimation 
of  the  remainder  of  battery  lifetime.  On  the  other  hand,  the  bat¬ 
tery  aging  also  induces  the  loss  of  useful  capacity  and  increase  of 
internal  resistance.  The  rated  capacity  measurement  must  apply  a 
full  charge-discharge  process  that  is  time-intensive  and  only  suit¬ 
able  for  testing  in  a  laboratory.  Instead,  the  battery  impedance 
can  be  either  directly  analyzed  by  a  spectroscopic  instrument  or 
be  estimated  through  the  transient  response.  For  instance,  some 
prevailing  methods  have  analyzed  the  battery  voltage  profile  by 
detecting  the  drop  depth  in  a  discharge  operation  [14,15]  or  sub¬ 
jecting  the  battery  to  a  specific  load  demand,  such  as  that  seen 
during  engine  cranking  [16,17],  in  order  to  determine  the  SoH. 
Prior  studies  [9,18-20]  have  also  demonstrated  that  the  accuracy 
of  the  SoH  determination  can  be  enhanced  though  a  more  elabo¬ 
rated  impedance  analysis.  By  employing  battery  electrical  models, 
the  internal  resistance  regarded  as  the  low  frequency  impedance 
can  be  estimated  by  exploiting  a  sliding-mode  control  technique 
[21],  or  an  extended  Kalman  filter  [12]. 

In  fact,  the  applicability  of  the  existing  SoC  or  SoH  determining 
methods  could  be  confined  to  some  specific  types  of  batteries.  For 
example,  the  internal  resistance  of  lithium-ion  battery  is  insensi¬ 
tive  to  SoC  variation  within  the  range  that  is  safe  for  use.  Yet,  the 
Coulomb  counting  or  OCV  measurements  have  no  such  limitation. 
Technically,  the  OCV  is  generally  measured  when  the  internal  ther¬ 
mal  stability  of  the  battery  has  been  reached.  To  do  this,  the  battery 
should  be  discharged  to  each  SoC  set-point  and  then  rested  for  a 
long  time  prior  to  measurement.  To  enhance  the  efficiency  of  the 
OCV  measurement,  a  statistical  analysis  based  on  the  data  obtained 
from  a  fast  discharge-charge  process  in  laboratory  has  been  pre¬ 
sented  [22].  In  practice,  a  complex  spectroscopic  analysis  of  battery 
impedance  is  unnecessary  as  long  as  the  internal  resistance  can  be 
estimated  instead.  Furthermore,  a  prior  study  has  illustrated  how 
to  determine  the  capacity  fade  of  battery  by  identifying  parameter¬ 
ized  curves  of  the  internal  resistance  with  respect  to  the  different 
working  cycles  based  on  a  simple  equivalent  circuit  model  (ECM) 
[19].  In  summation,  the  OCV  and  internal  resistance  parameters 
allow  for  determination  of  the  SoC  and  SoH,  respectively. 

In  this  study,  a  novel  method  for  online  estimating  of  lithium-ion 
battery  OCV  and  internal  resistance  simultaneously  is  presented. 
This  method  is  conceived  under  following  requirements  for  EV 
application:  (1)  online  estimation,  (2)  simple  mathematical  opera¬ 
tion,  (3)  no  need  of  a  specific  load  pattern  or  system  excitation,  (4) 
no  prior  time-intensive  training  process,  and  (5)  easy  implementa¬ 
tion.  To  this  end,  a  model-based  estimation  scheme  is  developed 
using  an  adaptive  control  approach  that  is  primarily  utilized  to 
design  controllers  for  nonlinear  or  linear  systems  with  uncertain¬ 
ties.  According  to  a  simple  ECM,  the  dynamics  of  a  lithium-ion 
battery  is  described  in  state-space  formulation,  and  the  estimated 
parameters  of  the  system  appear  linearly,  which  enables  utiliza¬ 
tion  of  adaptive  control  techniques.  The  persistent  excitation  (PE) 
required  for  assuring  parameter  convergence  in  an  adaptive  con¬ 
trol  system  is  shown  to  be  relaxed  so  that  the  convergent  time  that 
is  required  can  be  reduced.  The  proposed  estimation  system  is  the 
enhancement  of  the  previous  work  [23],  which  is  concerned  with 
a  method  of  estimation  for  battery  parameters. 


2.  Mathematical  battery  model 

Several  existing  battery  models  for  characterizing  battery 
electrochemical  behavior  have  been  constructed  by  using  math¬ 
ematical  approaches  from  either  macroscopic  or  microscopic 
perspectives.  For  microscopic  models,  partial  differential  equations 
are  entailed  for  modeling  the  electrochemical  interaction  between 
two  electrodes  and  electrolytes  in  order  to  determine  the  battery’s 
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Fig.  1.  A  generalized  ECM  for  lithium  batteries. 


effective  capacity,  current-voltage  relationship,  and  heat  gener¬ 
ation  [24,25].  However,  such  mathematical  models  are  generally 
complex  and  need  a  numerical  method  to  solve  the  problems, 
which  are  subject  to  initial  and  boundary  conditions.  On  the  other 
hand,  the  electrical  models  that  regard  the  battery  cell  as  a  lumped 
system  have  been  widely  introduced  to  capture  the  dynamic  char¬ 
acteristics  in  terms  of  current  and  voltage  in  order  to  ease  the 
computation  [26,27].  These  have  been  demonstrated  to  agree  with 
the  dynamical  response  in  discharge  or  charge  operation  within  the 
order  of  5%  error,  which  is  acceptable  in  terms  of  an  engineering 
viewpoint.  Due  to  the  lower  amount  of  required  computation  rel¬ 
ative  to  the  microscopic  models,  they  are  capable  of  determining 
the  battery  status  online  via  the  designed  algorithm. 

The  electrical  model  described  by  a  circuit  that  comprises  the 
basic  elements,  such  as  resistor,  inductor,  and  capacitor,  is  called  an 
ECM.  Those  components  are  configured  to  properly  match  the  mea¬ 
sured  impedance  or  transient  response  via  system  identification. 
Typically,  the  architecture  of  the  circuit  is  composed  of  a  funda¬ 
mental  ohmic  resistor  and  one  or  more  RC  networks  connected  in 
series  to  simulate  both  the  transient  and  steady  responses  of  the 
battery.  In  order  to  capture  the  nonlinear  characteristics  of  the  bat¬ 
tery,  the  resistance  of  each  identified  resistor  is  usually  a  nonlinear 
or  piece-wise  linear  function  with  respect  to  SoC  and  tempera¬ 
ture,  depending  on  the  types  of  batteries  used.  One  of  the  ECMs 
used  to  simulate  cell  performance  is  illustrated  in  Fig.  1,  where 
an  ohmic  resistor  with  resistance  RSl  an  RC  network  (Rt//Q),  and  a 
DC  source  with  voltage  voc  (OCV)  that  is  function  of  SoC  are  con¬ 
nected  in  series.  vb  is  defined  as  the  battery  terminal  voltage  and 
ib  is  the  outflow  current.  It  should  be  noted  that,  for  simplicity,  the 
above  parameters  are  generally  assumed  to  be  independent  of  cur¬ 
rent  direction.  In  other  words,  there  is  no  hysteretic  behavior  being 
considered  in  Fig.  1.  Accordingly,  the  impedance  between  the  DC 
source  and  positive  terminal  can  be  given  as 


z(jco)  =  Rs  + 


Rt 

jCtRtco  + 1 


(1) 


where  co  is  the  frequency  of  ib  or  vb.  The  parameters  in  Eq.  (1)  can 
be  identified  through  spectroscopic  analysis,  which  is  typically  an 
approach  applied  in  a  laboratory.  Alternatively,  it  is  more  straight¬ 
forward  to  use  battery  voltage,  current,  and  temperature  such  that 
an  online  estimation  can  be  realized.  To  this  end,  the  ECM  is  consid¬ 
ered  here  for  the  development  of  an  estimation  algorithm.  It  is  not 
only  applicable  to  the  scale  of  a  cell,  but  also  to  that  of  a  pack.  For 
instance,  a  more  precise  ECM  version  with  one  resistor  and  three 
RC  networks  in  series  was  utilized  for  modeling  a  battery  pack  to 
enhance  battery  modeling  [28]. 

Using  Kirchhoff  s  law,  the  dynamics  of  the  ECM  shown  in  Fig.  1 
can  be  expressed  as 


^=-ckVc+hib 


(2) 


Vb  =  Voc  -  Rsib  -  Vc 


(3) 


where  vc  is  defined  as  the  voltage  across  the  RC  network,  as  seen 
in  Fig.  1,  and  cannot  be  determined  because  the  voltage  relating  to 
Rs  is  unknown.  It  can  be  seen  that  Eqs.  (2)  and  (3)  are  the  state- 
space  formulation,  where  vc  can  be  regarded  as  an  estimated  state 
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for  an  observer  or  Kalman  filter  design  [12]  as  long  as  Rs,  Rt,  Ct,  and 
voc  are  previously  known.  However,  this  is  not  the  case  because  all 
the  above  parameters  along  with  voc  are  still  uncertain.  In  addi¬ 
tion,  more  precisely,  they  should  be  variant  with  respect  to  SoC, 
temperature  (T),  and  usage  history  (h),  i.e.: 

p  =  p(SoC,  T.h),  pT=[Rs  Rt  Ct  Foe]  (4) 

where  p  is  a  parametric  vector  with  the  inclusion  Rs,  Rtl  Ct,  and  voc. 
All  of  the  battery  parameters  in  (4)  vary  slightly  within  90-10%  of 
the  SoC,  implying  that  they  can  be  estimated  by  any  methodology 
that  can  account  for  such  slowly  varying  behavior. 

Because  vc  is  neither  a  measurable  variable  nor  an  estimated 
parameter  following  the  above  formulation,  the  removal  of  the  vc 
term  from  Eqs.  (2)  and  (3)  is  desired  in  order  to  facilitate  the  algo¬ 
rithm  derivation.  To  this  end,  one  can  differentiate  vb  in  Eq.  (2)  to 
yield: 

vb  =  Voc  -  Rsib  -  Rsib  -Vc  =  &{Voc)x  -  v(Rs)xib  -  Rsib  ~  ”c  (5) 
where  the  differential  operator  cr(A)  and  vector  x  are  defined  as 


dp 
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Here,  p  is  any  parameter  in  p.  It  can  be  seen  that  the  first  two 
terms  on  the  right-hand  side  of  Eq.  (5)  are  generated  by  deriving 
the  higher-order  dynamics  of  the  battery  in  terms  of  the  output 
voltage  vb.  The  complicated  variation  of  a  parameter  is  interpreted 
as  multiplication  of  the  change  of  the  parameter  with  a  state  (SoC, 
T,  or  h),  and  the  change  of  the  state  with  time.  In  addition,  Eq.  (5) 
can  be  simplified  with  the  following  assumptions: 

(Al)  Small  battery  energy  is  consumed  or  regained  relative  to 
totally  useful  capacity.  For  example,  a  40  Ah  lithium-ion  cell, 
the  basic  energy  requirement  for  a  city-EV,  is  assumed  to 
conduct  1  C  (40  A)  discharge  (generally  less  than  this  value 
at  most  time).  The  variation  of  SoC  with  respect  to  time  is 
3SoC/9 1=-  40/(40  x  3600)  =  -0.00028,  which  is  very  small, 
i.e.,  dSoC/dt^O. 

(A2)  The  cell  temperature  must  be  monitored  and  controlled  at  a 
predetermined  level  by  a  BMS  in  order  to  avoid  a  thermal  run¬ 
away  problem  and  keep  thermal  stress  as  low  as  possible  for 
the  longevity  of  the  battery.  Relying  on  the  proper  design  of 
a  cooling  system/heater,  the  cell  temperature  rise/decrease 
should  be  slow.  Hence,  dT/dt^O  holds  for  normal  operating 
conditions. 

(A3)  Note  that  h  represents  a  long  term  usage  history  such  as  the 
counts  of  working  cycles  or  abuse.  Thus,  dh/dt^O  definitely 
holds. 


Accordingly,  it  follows  that  x  ~  0,  thereby  Eq.  (5)  can  be  rewrit¬ 
ten  as 

•  .  •  1  1  .  •  1  . 

Vb  ^  —RSlb  -VC  =  ~RSlb  +  Tttt  VC  ~  -pr  lb  =  ~RSlb  ~  7^ lb 

1  LfKf  L  t  L  t 

+  r  p  (^oc  —  Rsib  —  vb)  (6) 

Cf/vf 

_  j_  n  •  Rt  +  Rs  •  .  Voc 

Ct  Rt  b  s  b  CtRt  b  +  CtRt 


the  system  inputs  multiplied  by  the  coefficients  l/CtRt,  Rs,  and 
( Rt  +  Rs)ICtRt ,  respectively,  to  give  new  parameters  to  be  deter¬ 
mined.  Because  the  parameters  in  Eq.  (6)  appear  linearly,  it  is 
possible  to  employ  an  adaptive  control  approach  that  continuously 
updates  uncertain  parameters  before  the  tracking  error  approaches 
zero.  Recall  that  OCV  is  slowly  varying  within  the  battery’s  use- 
able  capacity.  The  region  between  10%  and  90%  of  SoC  is  generally 
unused  in  order  to  avoid  the  abusive  conditions  such  as  overcharge 
and  overdischarge.  Hence,  OCV  is  considered  as  one  of  the  param¬ 
eters  such  that  Eq.  (6)  can  be  written  as 

h  =  0TX  (7) 


where 


oT=[Ot  e2  e3  e4f 


p  Rt  +  Rs  1  Fqc 

.  s  ~QRr  QRt  CtRt  J  ’ 


XT  =  [  -ib  -ib  -vb  1  ]  . 

Note  that  voc  and  1 1  QRt  are  coupled  together  to  become  one  new 
parameter.  The  estimated  parameter  vector  in  Eq.  (7)  is  defined  as 


oT  =[9t  e2  03  e4f  = 

The  corresponding  predicted  state  vb  is  given  as 
vb  =  0  -\-  u 


Rt  +  Rs 

CcRt 


1  Foe 

QRt  QRt 


(8) 


(9) 


where 

XT  =  [  -ib  -ib  -vb  1  ] 

and  u  is  input  designed  for  regulation. 

The  design  of  an  adaptive  control  system  is  to  let  a  controlled 
state  follow  a  reference  signal.  In  our  case,  we  expect  that  the  track¬ 
ing  error  of  the  predicted  state  can  asymptotically  tend  towards 
zero,  i.e.: 


lime  =  lim (vb(t)  -  vb(t ))  =  0.  (10) 

t->  OO  t->0 o 

It  is  known  that  Eq.  (10)  is  also  necessary  for  the  convergence  of 
the  estimated  parameters.  By  utilizing  the  technique  used  for  the 
design  of  the  adaptive  controller  [29],  an  adaptation  law  of  the 
parameters  can  be  developed. 


Theorem.  For  a  battery  formulated  in  Eq.  (7),  its  estimated 
dynamics  is  given  by  Eq.  (9),  and  the  condition  expressed  by  Eq. 
(10)  can  be  achieved  by  the  adaptation  law: 


-gihe 

^2 

-g2he 

03 

-g3%e 

A. 

g4e 

together  with  the  designed  input: 


u  =  Xe 


(12) 


where  A  is  positive  semi-definite,  i.e.: 


by  substituting  Eqs.  (2)  and  (3)  into  Eq.  (5).  Obviously,  the  term  vc  A  >  0, 
has  been  removed  and  a  new  variable  ib  obtained.  Yet,  ib  can  be 

easily  computed  by  the  differential  of  measured  ib.  and  the  adaptation  gains  are  positive  definite: 


3.  Main  results 


gi  >  0,  i  =  1,2,  3, 4. 


3.1.  Adaptive  algorithm 


Proof.  Invoking  the  Lyapunov-based  design  procedure,  consider 
the  Lyapunov  function: 


According  to  Eq.  (6),  the  dynamics  of  the  battery  is  reinterpreted 
in  terms  of  vb.  On  the  other  hand,  voc,  ib,  and  ib  are  regarded  as 


V(.vb,d)='-e2+'-tfJ2~° 


(13) 
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where  E  is  a  positive  definite  matrix  with  proper  dimensions.  The 
estimated  error  vector  0  defined  as  the  difference  between  actual 
and  estimated  parameters  is 

oT  =  oT-dT=  [g,  e2  03  04].  (i4) 

Note  that  V  is  a  positive-definite  function  constituted  by 
a  quadratic  polynomial  in  e  and  6.  Hence,  if  V  monotonically 
decreases,  Eq.  (10)  is  achieved.  It  is  sufficed  that  its  derivative  shall 
be  negative  semidefinite.  To  this  end,  differentiating  Eq.  (13)  along 
the  trajectory  of  the  system,  it  follows  that 

^  _  dV_dVb  dVdO 
-  dh~dF  +  ag  dt 
=  ee  +  6  ^0 
=  {vb-vb)e  +  {0T  -0 
=  [0TX  -  eTX  -  u)e  +  [0T  -  0  )£§ 

=  -0i'ib-02ib-(03+03)(e  +  vb)  +  03Vb  +  04-u  e 

.  -r  .  J  (15) 

+(0T-0)E 0 

=  [-Wh  ~  hk  ~  03%  +  04  -  03e  -  u]  e  +  (6»T  -  0 
=  {XTe-0  e  +  u)e  +  0TJ20 

=  -(03  +  Ve2+0TYl0  ^ 

<0 

provided  that 
u  =  Ae,  A  >  0, 


V 

II 

<X 

Tw 

II 

0  = 

02 

b3 

-giibe 

-g3%e 

A. 

g4? 

by  letting  ^=diag((l/g,)>0,  i  =  1,2, 3),  and 

0  =  a(0)x  ^  0  (16) 

is  used  in  accordance  with  (A1)-(A3).  Here  r t  =  QRt  is  the  polar¬ 
ization  time-constant  with  regard  to  battery  relaxation  behavior. 
Particularly,  the  designed  gain  A  can  be  selected  to  be  larger  in 
order  to  dominate  the  convergent  rate  so  as  to  shorten  estimation 
time.  By  applying  the  LaSalle-Yoshizawa  theorem  [30]  to  Eq.  (15), 
Eq.  (10)  follows.D 

3.2.  Convergence  of  parameters 


It  can  be  found  from  Eq.  (9)  that  6\ ,  when  coupled  with  ib,  makes 
itself  more  sensitive  to  high  dynamic  operation,  e.g.,  EV  traction 
as  well  as  regeneration,  in  comparison  with  the  other  parameters, 
implying  that  the  convergent  rate  of  will  be  much  faster  than 
those  of  the  other  parameters.  As  a  consequence,  the  condition: 

0,^0  (17) 


will  be  easily  achieved  in  a  short  time.  This  phenomenon  can  be 
observed  in  tentative  simulations  in  the  following  section. 
According  to  Eqs.  (10)  and  (17),  we  have 

lim  (vb(t)  -  %(t))  =  lim  (0TX(t)  -  -  u(t)) 

t->  OO  t — >  OG 


=  -ib{t)  lim 0i (t)  -  ib{t) lim 02(t) 

t — >  oo  t—^oo 

-vb{t)  lim  03(t)  +  lim94(t)-(03(t) 

t — ^  oo  t—>  oo 

+A)lime(t) 

t — ^  oo 

=  -ib(t) lim 02(t)  -  Vb(t)lim03(t) 

t — >  oo  t — ^  OO 

+  limS4(t). 


(18) 


Eq.  (18)  should  equal  to  zero  when  the  stability  condition  has  been 
achieved,  i.e.: 

-ifc(t) lim 02(O  -  vb{t)\im03{t)  +  lim04(t)  =  0.  (19) 

t — >■  oo  t — ^  OO  t->0 o 


On  the  other  hand,  let  Eq.  (3)  be  multiplied  by  1  ICtRt  on  both 
sides  such  that 


Pb  Poc  Rs  •  1  ^  Poc  Rs  +  Rt  ■ 

RtCt  =  RtQ  ~  lCctlb  ~  lCctVc  ~^cTlb 


(20) 


Consider  a  situation  where  the  RC  network  in  the  ECM  has  a  low 
time-constant  (rt  =  RtQ)  such  that  it  acts  like  a  pure  resistor,  i.e.: 


Pc^Rtib  (21) 

so  as  to  make  it  possible  to  rewrite  Eq.  (20)  as 

-02ib  ~  03Pb  +  0A  *  0.  (22) 

Multiplying  Eq.  (22)  by  the  constant  a  -  1  and  then  adding  it  to  Eq. 
(18)  yields: 

-ib{t)\im(a02  -  02(t))  -  vb{t)\im(a03  -  03(t)) 

t— ►  OO  t — >  OG 


+  lim(a04-04(t))«O.  (23) 

t — ^OO 


It  is  illustrated  that  the  adaptation  algorithm  achieves  Eq.  (10) 
and  forces  all  estimated  parameters  to  approach  the  values  in  pro¬ 
portion  to  their  actual  values  scaled  by  an  unknown  constant  a , 
i.e.: 


Once  0  converges  to  its  actual  value,  each  parameter  can  be 
evaluated  via  Eq.  (8).  It  should  be  noted  that  Eq.  (15)  assures 
asymptotical  stability  of  the  state  error  e,  rather  than  the  esti¬ 
mated  parametric  error  0,  as  it  is  purposely  cancelled  to  achieve 
sufficient  system  stability.  As  a  general  shortcoming  of  an  adaptive 
control  system,  the  PE  is  typically  needed  to  assure  the  parameter 
stability  where  the  PE  is  a  type  of  design  input  needed  to  generate 
rich  enough  data  for  system  identification  [31  ].  Because  the  battery 
response  is  primarily  determined  by  current  command  of  the  trac¬ 
tion  motor,  the  driving  behavior  may  not  produce  a  suitable  current 
pattern  to  satisfy  the  PE.  As  a  result,  the  estimated  parameters  may 
either  take  more  time  for  convergence  or  become  trapped  at  an 
incorrect  equilibrium  point.  To  evaluate  the  above  problems,  the 
stability  condition  of  state  error  e  should  be  analyzed  as  it  appears 
in  each  term  of  the  adaptation  law  (see  Eq.  ( 1 1 )),  i.e.,  the  parameters 
stop  updating  as  long  as  the  state  error  approaches  zero. 


lim  Gi(t)  =  aGi,  i  =  2,3,4,  (24) 

t— >oo 


which  is  caused  somehow  by  initial  condition,  driving  pattern,  and 
adaptation  gain  design.  Unfortunately,  this  is  the  case  in  both  sim¬ 
ulations  and  experiments,  as  will  be  illustrated  in  a  later  section. 
Although  actual  values  of  the  model  parameter  set  are  difficult  to 
obtain  in  a  straightforward  manner,  it  is  possible  to  evaluate  all 
battery  parameters  other  than  the  capacitance  Q.  To  see  this,  by 
using  Eqs.  (17),  (20)  and  (24)  we  have 


Rs(t) 

RtU) 

t >oc(t) 


01  (t) 


hit) 

Ht) 


0i  (t) 


04  (0 

-  03(t)  - 


(25) 
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Fig.  2.  Scheme  of  SoC  and  SoH  determination  based  on  adaptive  control  approach. 


where  [rs  Rt  ] T  and  voc  can  be  used  to  determine  SoH  and  SoC, 
respectively. 

3.3.  Implementation  of  SoC  and  SoH  determination 

The  proposed  estimation  method  provides  online  monitoring 
of  the  battery’s  internal  resistance  and  OCV  for  the  determination 
of  SoH  and  SoC,  respectively.  During  the  estimation  procedure, 
a  front-end  process  acquires  current,  voltage,  and  temperature 
via  direct  measurements.  These  signals  or  data  may  contain  mea¬ 
suring  noise  such  that  a  filtering  process  might  be  necessary  to 
improve  the  estimation  reliability.  Later,  the  filtered  data  are  sent 
to  the  developed  estimation  process  for  extracting  battery  internal 
resistance  and  OCV.  The  fact  that  the  internal  resistance  gradually 
increases  as  the  battery  decays  can  be  used  to  convert  the  estimated 
resistance  data  to  a  human-sense  unit  such  as  battery  age  or  healthy 
ratio.  On  the  other  hand,  the  OCV  curve,  which  decreases  propor¬ 
tionately  with  the  battery’s  usable  energy,  is  used  to  determine  the 
SoC.  In  addition,  the  battery’s  working  temperature  effects  on  the 
electrical  characteristics  shall  be  considered  in  the  above  SoC  and 
SoH  conversions.  In  order  to  reduce  the  estimation  susceptibility 
to  the  measuring  precision  or  disturbance,  a  filter  with  high  band¬ 
width  is  introduced  to  process  the  estimated  parameters  prior  to 
the  conversion.  A  detailed  scheme  for  the  SoC  and  SoH  determina¬ 
tion  is  illustrated  in  Fig.  2.  The  relevant  processes  are  summarized 
as  follows. 

1.  Measurement:  A  BMS  typically  provides  cell  voltage  and  temper¬ 
ature  monitoring  by  integrating  dedicated  measuring  modules. 
Battery  terminal  voltage  can  be  obtained  by  summing  the  cell 
voltages  in  series.  A  Hall-effect  current  sensor  is  deployed  in  the 
series  loop  for  current  management.  Relying  on  the  general  BMS, 
the  data  requested  for  the  estimation  are  acquired. 

2.  Filter  I  and  filter  2:  The  bandwidth  of  the  filter  1  shall  be  properly 
chosen  to  filter  out  only  measurement  noise  without  eliminat¬ 


ing  battery  dynamics.  However,  the  time  constant  of  the  filter  2 
should  be  selected  to  be  sufficiently  large  to  enhance  the  robust¬ 
ness  of  the  estimation  against  measuring  disturbance.  Also,  the 
parameters  of  the  filters  are  initially  configured  and  held  invari¬ 
ant. 

3.  Battery  model:  The  ECM  given  by  Eqs.  (2)  and  (3)  is  reformulated 
into  Eq.  (7).  It  is  assumed  that  Eq.  (7)  can  catch  dynamics  of  bat¬ 
teries  used  for  EV  applications  and  can  be  observed  by  Eq.  (9) 
along  with  the  estimated  parameter  vector  Eq.  (8). 

4.  Adaptive  algorithm:  The  adaptation  law  given  by  Eqs.  (11)  and 
(12)  is  processed  to  update  the  estimated  parameters. 

5.  Parameter  extraction:  The  internal  resistance  and  OCV  are 
extracted  via  Eq.  (25). 

6.  SOC  and  SOH  determination:  The  lookup  tables  define  data  map¬ 
pings  between  battery  parameter  and  working  status  models  in 
consideration  with  thermal  effect.  They  are  made  through  long¬ 
term  tests  in  a  laboratory  or  possibly  provided  by  battery  makers 
for  such  application. 

4.  Simulation  and  experimental  results 

The  simulations  and  experiments  are  performed  to  validate  the 
effectiveness  of  the  proposed  method  on  the  estimation  of  bat¬ 
tery  internal  resistance  and  OCV.  Two  simulations  established  in 
Simulink  are  used  to  verify  the  capability  of  determining  the  param¬ 
eters  in  ECM  accurately.  The  OCV  is  set  to  be  fixed  and  slowly 
varying  with  respect  to  SoC  in  Simulations  1  and  2,  respectively. 
Furthermore,  two  EV  battery  packs  constructed  from  different 
types  of  lithium  cells  for  a  light  EV  [32]  are  utilized  in  two  exper¬ 
iments  to  justify  the  applicability.  To  mimic  a  scenario  in  which 
the  battery  pack  supplies  motive  power  to  the  light  EV  running 
in  an  urban  area,  the  American  FTP-75  driving  cycles  are  adopted 
to  generate  the  demanded  power  pattern  by  means  of  ADVISOR,  a 
vehicle  simulator.  Automated  cell-test  equipment  commands  the 
attached  battery  pack  to  discharge  or  charge  following  a  power 


Table  1 

Battery  parameters  and  adaptation  gains. 


Group 

Symbol 

Simulation  1 

Simulation  2 

Experiment  1 

Experiment  2 

Battery 

Voc  (V) 

12.6 

12. 6a 

70.7a 

69.3a 

parameters 

Rs  (m£2) 

0.1 

0.1 

76.82b 

10.22b 

Rt  (m£2) 

0.05 

0.05 

- 

- 

Q  (F) 

5 

5 

- 

- 

Adaptation 

X 

0.2 

0.2 

0.4 

0.4 

gains 

Si 

0.001 

0.001 

0.001 

0.0004 

&2 

0.005 

0.005 

0.5 

0.5 

S3 

0.05 

0.05 

0.8 

0.8 

§4 

1 

1 

1 

1 

a  The  OCV  is  varying  with  the  change  of  the  capacity.  The  value  shown  here  represents  its  initial  state  in  simulation  or  experiments. 
b  The  ohmic  resistance  is  measured  by  an  impedance  measuring  instrument  as  its  actual  value. 
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Time  /  s 


Fig.  3.  Measurements  and  estimated  state  error  in  Simulation  1 :  (a)  terminal  battery 
voltage,  (b)  battery  outflow  current,  and  (c)  estimated  state  error. 

pattern  in  a  laboratory.  The  voltage  and  current  of  the  battery  pack 
are  measured  via  the  equipment  and  recorded  in  its  remote  PC. 
Because  the  scheme  works  in  an  open-loop  manner,  the  estima¬ 
tion  process  can  be  performed  offline  without  loss  of  reality.  To 
this  end,  the  recorded  data  are  loaded  into  the  Simulink  environ¬ 
ment  to  simulate  the  real  estimation  condition.  In  the  following 
simulations  and  experiments,  the  measured  states  are  battery 


working  voltage  vb  and  current  ib.  The  estimated  battery  volt¬ 
age  vb  is  produced  by  the  ECM  with  estimated  model  parameters 
i  =  1~4  as  well  as  the  measured  current.  The  primary  estimated 
state  for  SOC  is  OCV  voc,  and  the  primary  estimated  parame¬ 
ters  for  SOH  are  internal  resistance  Rs  and  Rt.  voc,  Rs,  and  Rt  can 
be  extracted  from  the  model  parameters.  The  system  parame¬ 
ters  required  for  the  simulations  and  experiments  are  specified  in 
Table  1. 


4.1.  Simulation  1 

The  objective  of  this  simulation  is  to  verify  the  accuracy  of 
the  proposed  method  when  estimating  the  internal  resistance 
and  OCV.  A  battery  exactly  described  by  Eqs.  (2)  and  (3)  is 
assumed  to  discharge  and  charge  subject  to  a  random  load  as 
depicted  in  Fig.  3(a).  The  load  voltage  is  equal  to  the  battery 
terminal  voltage  and  determines  the  current  profile  as  shown 
in  Fig.  3(b),  where  the  positive  value  denotes  outflow  current. 
The  estimated  parameters  are  initially  set  at  0i(O)  =  fls(0)  =  0.2  £2, 
e2{ 0)  =  (Rt  +  R07ctRt(0)  =  o.oi  F"1 . 03{0)  =  f/QRt(0)  =  6  s”1 .  and 
04(0)  =  t'Oc/QRt(0)  =  40  Vs-1.  The  trajectory  of  the  estimated  error 
is  depicted  in  Fig.  3(c)  and  converges  to  zero  within  10  s.  The  esti¬ 
mated  parameters  and  their  actual  values  are  depicted  in  Fig.  4, 
where  the  blue  and  black  lines  denote  the  actual  and  estimated  val¬ 
ues,  respectively.  It  can  be  found  that  only  6\  agrees  with  the  actual 
value  of  0.1  £2,  reflecting  the  observation  made  in  Eq.  (17),  and  the 
others  stop  updating  after  the  convergence  of  the  estimated  error. 
The  resulting  internal  resistance  and  OCV  are  depicted  in  Fig.  5  and 
both  tend  to  their  actual  values  within  40  s  despite  the  accuracy  and 
traces  of  the  estimated  parameters.  This  demonstrates  the  effec¬ 
tiveness  of  the  proposed  method  when  applied  to  the  ECM  along 
with  the  invariant  parameters. 


Fig.  4.  Trajectories  of  estimated  parameters  in  Simulation  1:  (a)  9\,  (b)  02,  (c)  O3,  and  (d)  #4;  the  blue  and  black  lines  denote  actual  and  estimated  values,  respectively.  (For 
interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 
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Fig.  5.  Estimated  internal  resistance  and  OCV  in  Simulation  1 :  (a)  Rs,  (b)  voc,  and  (c) 
Rt- 
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Fig.  6.  Measurements  and  estimated  state  error  in  Simulation  2 :  (a)  terminal  battery 
voltage,  (b)  battery  outflow  current,  and  (c)  estimated  state  error. 


4.2.  Simulation  2 

A  ECM  where  the  OCV  is  monotonically  varying  throughout  the 
entire  operation  is  considered  in  this  simulation.  Other  than  the 
OCV,  the  parameters  and  their  initial  values  are  invariant,  simi¬ 
lar  to  the  ones  set  in  Simulation  1.  A  random  load  profile  used  to 
excite  the  system  is  also  the  same  and  is  depicted  in  Fig.  6(a).  Due 


to  the  rise  of  the  OCV,  the  battery  supplies  more  and  regains  less 
current  during  the  discharge  and  charge  operations,  respectively, 
as  depicted  in  Fig.  6(b)  when  comparison  with  Fig.  3(b).  It  can  be 
observed  from  Fig.  6(c)  that  the  estimated  error  converges  to  zero 
within  10  s.  The  estimated  parameters  are  depicted  individually  in 
Fig.  7  for  comparison  to  their  actual  values,  and  the  following  brief 
conclusions  are  drawn. 


Fig.  7.  Trajectories  of  estimated  parameters  in  Simulation  2:  (a)  0\ ,  (b)  02,  (c)  O3,  and  (d)  #4;  the  blue  and  black  lines  denote  actual  and  estimated  values,  respectively.  (For 
interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 
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Fig.  8.  Estimated  internal  resistance  and  OCV  in  Simulation  2:  (a)  Rs,  (b)  voc,  and  (c) 
Rt- 


Fig.  9.  Measurements  and  estimated  state  error  in  Experiment  1 :  (a)  terminal  bat¬ 
tery  voltage,  (b)  battery  outflow  current,  and  (c)  estimated  state  error. 


(i)  The  estimated  ohmic  resistance  identical  to  6\  converges  to  its 
actual  value  as  fast  as  the  case  in  Simulation  1,  as  shown  in 
Fig.  7(a).  It  is  evident  that  the  condition  set  forth  in  Eq.  (17)  is 
true  in  either  two  cases. 

(ii)  Comparatively,  the  other  estimated  parameters,  diverge 
from  the  actual  profiles  as  illustrated  in  Fig.  7(b)-(d).  How¬ 
ever,  the  resulting  internal  resistance  and  OCV  still  tend  to 
their  actual  values  after  10  s  as  shown  in  Fig.  8(a)-(c).  It  can 


be  inferred  from  Figs.  7  and  8  that  the  multiplier  a  is  varying 
during  the  entire  process,  rather  than  being  held  constant  as  in 
Simulation  1. 

Simulations  1  and  2  verify  the  proposed  method  based  on  the 
assumptions  (A1)-(A3)  and  can  accurately  estimate  the  internal 
resistance  and  OCV  of  the  battery  that  can  be  exactly  modeled  by 
Eqs.  (2)  and  (3). 


Fig.  10.  Trajectories  of  estimated  parameters  in  Experiment  1 :  (a)  9\,  (b)  O2,  (c)  O3,  and  (d)  64. 
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Fig.  11.  Trajectories  of  estimated  internal  resistance  and  OCV  in  Experiment  1:  (a) 
RSl  (b)  Voc,  and  (c)  Rt. 

4.3.  Experiment  t 

Practically,  the  ECM  employed  here  is  simplified  to  reduce  the 
order  of  model  complexity  and  facilitate  the  development  of  the 
estimation  algorithm.  The  accuracy  affected  by  unmodeled  dynam¬ 
ics  should  be  investigated  by  using  a  real  battery.  To  this  end, 
lithium  iron  phosphate  cells  are  connected  in  series  to  reach  a  nom¬ 
inal  voltage  of  70  V  to  fulfill  the  voltage  requirement  of  a  light  EV. 
An  experiment  is  established  to  let  the  automated  cell-test  machine 
manage  the  discharge  and  charge  operation  in  accordance  with 
the  converted  power  pattern.  Meanwhile,  the  ambient  tempera¬ 
ture  of  the  battery  pack  is  managed  by  an  air  conditioner  to  remove 
the  accumulated  heat  within  the  cells  during  operating.  Thus,  the 
cells’  temperatures  would  appear  in  slowly  varying  manner  so  as 
to  meet  the  assumption  (A2).  Throughout  one  cycle  of  the  FTP-75, 
the  battery  voltage  and  outflow  current  are  depicted  in  Fig.  9(a) 
and  (b),  respectively.  The  voltage  fluctuates  between  73  and  64  V 
and  the  maximum  discharge/charge  current  varies  up  to  95/40  A. 
The  initial  values  for  the  estimated  parameters  are  0i(O)  =  0.1  £2, 
@2(0)  =  20 F-\  03(O)  =  100 s-1,  and  @4(0)  =  7067.2  Vs”1.  As  illus- 
trated  in  Fig.  9(c),  the  tracking  error  converges  to  zero  after  5  min. 
With  proper  gain  selection  and  prior  knowledge  of  initial  paramet¬ 
ric  values,  faster  convergence  can  be  achieved.  The  trajectories  of 
the  estimated  parameters  are  depicted  in  Fig.  10(a)-(d).  As  shown 
in  Fig.  10(a)-(c),  the  parameters  0\,02,  and  03  converge  to  steady 
values  of 0.063, 7.5,  and  1 00.5,  respectively,  with  small  fluctuations. 
Unlike  the  others,  the  parameter  04  follows  a  descending  manifold, 
as  illustrated  in  Fig.  10(d),  as  the  OCV  drops  as  a  function  of  the 
decrease  of  battery  stored  energy. 

The  resulting  internal  resistance  and  OCV  are  depicted  in  Fig.  11. 
The  ohmic  resistance  trajectory  equals  to  the  parameter  as 
shown  in  Fig.  11(a).  The  polarized  resistance  tends  to  0.01  £2  and 
the  OCV  descends  slowly  from  its  initial  value  of  70.75  to  70.15  V 
at  the  end  of  the  process.  In  order  to  enhance  the  reliability  of  SoC 
and  SoH  determination,  the  estimated  internal  resistance  ( Rs  +  Rt) 
and  OCV  are  processed  through  the  filter  2  to  become  smooth,  as 
shown  in  Fig.  1 2.  Rs  +  Rt,  regarded  as  low  frequency  impedance  that 
is  approximated  as  75  m£2,  is  close  to  76.82  m£2,  the  value  obtained 
via  an  impedance  measuring  instrument,  as  shown  in  Table  1.  It  is 
shown  in  Fig.  12(b)  that  the  estimated  OCV  declines  slowly  in  the 


(b) 


Fig.  12.  Filtered  estimated  internal  resistance  and  OCV  in  Experiment  1 :  (a)  internal 
resistance  and  (b)  OCV. 


process  except  for  the  first  5  min  of  the  adaptation  process.  The 
accuracy  analysis  of  the  OCV  estimation  will  be  verified  in  the  next 
experiment. 

4.4.  Experiment  2 

A  lithium  polymer  battery  pack  is  constructed  to  reach  a  nomi¬ 
nal  voltage  of  68  V  and  a  higher  energy  density  than  the  one  used  in 
the  first  experiment  in  order  to  investigate  the  compatibility  of  the 
proposed  method  with  the  different  types  of  EV  batteries.  Again,  the 
battery  pack  is  managed  to  discharge  and  charge  in  accordance  with 
the  same  power  profile  and  thermal  environment  used  in  Experi¬ 
ment  1 .  The  estimated  parameters  are  initially  set  at  0\  (0)  =  0.05  £2, 
@2(0)  =  15 F— ^ ,  03(O)  =  100s-1,  and  @4(0)  =  6875 Vs"1.  The  bat- 
tery  working  voltage,  current,  and  estimated  error  are  depicted  in 
Fig.  13.  The  voltage  fluctuates  between  70  and  65  V  and  its  range  is 
larger  than  that  of  the  pervious  battery  pack  due  to  their  different 
electrochemical  properties.  The  estimated  error  has  been  attenu¬ 
ated  to  near  zero  after  450  s.  The  estimated  parameters  depicted  in 
Fig.  14  have  trajectories  similar  to  those  in  Simulation  2.  It  is  noted 
in  Fig.  14(a)  and  (b)  that  the  parameters  tripped  off  their  steady  val¬ 
ues  at  the  end  of  the  process.  As  shown  in  Fig.  13,  between  1 ,300  and 
1,350  s,  the  voltage  spikes  and  outflow  current  have  the  same  side, 
which  disagrees  with  the  normal  voltage-current  behavior  because 
the  internal  resistance  would  cause  excess  voltage  drop/bounce  on 
the  discharge/charge  operation.  This  fault  might  be  caused  by  a 
measurement  or  data  transmission  error  that  can  be  regarded  as  a 
disturbance.  The  resulting  ohmic  resistance,  OCV,  and  polarization 
resistance  are  plotted  in  Fig.  1 5(a)-(c),  respectively.  It  can  be  shown 
that  the  disturbance  affects  the  estimated  internal  resistance  more 
than  the  OCV.  Rs  +  Rt  and  estimated  OCV  processed  through  the  fil¬ 
ter  2  are  shown  in  Fig.  16.  After  450  s,  Rs  +  Rt  converges  to  a  steady 
value  of  20  m£2,  and  the  estimated  OCV  approaches  a  declined  man¬ 
ifold.  Also,  the  robustness  of  the  proposed  scheme  is  improved  by 
examining  the  Rs  +  Rt  curve  at  approximately  1,320  s.  In  actual¬ 
ity,  the  internal  resistance  measured  by  the  impedance  measuring 
instrument  is  10.22  m^,  as  shown  in  Table  1,  and  is  close  to  the 
estimated  ohmic  resistance  rather  than  low  frequency  impedance. 
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Fig.  13.  Measurements  and  estimated  state  error  in  Experiment  2:  (a)  terminal 
battery  voltage,  (b)  battery  outflow  current,  and  (c)  estimated  state  error. 

For  the  comparison  with  actual  OCV  at  the  beginning  of  estima¬ 
tion  process,  the  parameters  are  initially  set  to  be  equal  to  those 
obtained  from  the  previous  results,  i.e.,  @i(0)  =  0.01  £2,  6 2(0 )  = 
2F-1,03(O)  =  100s_1,  and  04(O)  =  6882.6  Vs-1.  By  reusing  the  bat- 


Fig.  15.  Trajectories  of  estimated  internal  resistance  and  OCV  in  Experiment  2:  (a) 
Rs,  (b)  voc,  and  (c)  Rt. 

tery  voltage  and  current  data,  the  estimated  OCV  with  respect  to 
depth  of  discharge  (DoD)  after  the  filtering  process  is  depicted  in 
Fig.  17.  On  the  other  hand,  the  actual  OCV  is  measured  after  a  pro¬ 
cess  in  which  the  battery  pack  is  discharged  with  a  low  current  of 


Fig.  14.  Estimated  parameters  in  Experiment  2:  (a)  0\,  (b)  O2,  (c)  O3,  and  (d)  64. 
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(i)  For  a  battery  exactly  modeled  by  the  ECM  shown  in  Fig.  1,  the 
proposed  method  can  be  performed  online  and  can  accurately 
estimate  the  internal  resistance  and  OCV. 

(ii)  For  an  EV  battery  pack  constructed  by  either  lithium  polymer 
or  lithium  iron  phosphate  type  of  cells,  the  proposed  scheme 
can  be  performed  online  and  can  reliably  estimate  the  ohmic 
resistance,  polarized  resistance,  and  OCV  with  a  small  error  so 
as  to  determine  the  SoC  and  SoFI. 

(iii)  The  proposed  scheme  is  robust  to  the  disturbance  caused  by 
measuring  or  data  transmission  error. 


Fig.  16.  Filtered  estimated  internal  resistance  and  OCV  in  Experiment  2:  (a)  internal 
resistance  and  (b)  OCV. 


Fig.  17.  Comparative  results  of  estimated  and  actual  OCV  with  respect  to  DOD. 


5  A  until  the  2%  DoD  is  met  and  then  rested  for  30  min.  This  pro¬ 
cess  is  repeated  four  times,  and  the  measured  voltage  is  plotted  in 
Fig.  1 7  by  a  solid  line.  The  actual  OCV  at  the  2%,  4%,  6%,  and  8%  DoD 
are  indicated  in  comparison  with  the  estimated  ones.  It  is  found  that 
the  maximum  estimated  error  is  ±0.1  V  in  OCV  or  ±0.5%  in  DoD  (or 
SoC).  The  overall  estimated  error  for  the  lithium  polymer  batteries 
without  considering  thermal  effect  is  about  ±1%  SoC  at  maxi¬ 
mum,  by  adding  maximum  estimation  discrepancy  (±0.1  V  OCV) 
at  the  same  SoC  (seeing  the  estimated  OCV  curve  at  3.3%  DoD  in 
Fig.  17). 

The  capability  of  the  proposed  scheme,  according  to  the 
simulation  and  experimental  results  are  summarized  as 
follows. 


5.  Conclusions 

In  this  study  we  proposed  a  technique  for  identifying  battery 
internal  resistance  and  OCV  regarding  SoFI  and  SoC  determination, 
respectively.  Based  on  an  ECM  simply  composed  of  one  resistor 
and  one  RC  network  connected  in  series  to  capture  the  essential 
dynamics  of  a  lithium-ion  battery,  an  online  estimation  method  of 
the  model  parameters  is  developed  by  using  an  adaptive  control 
approach  that  has  been  widely  applied  to  nonlinear  control  sys¬ 
tems  with  parameter  uncertainties.  The  adaptation  law  is  devised 
to  assure  that  estimated  state  error  asymptotically  decreases  to 
zero  by  invoking  the  Lyapunov  stability  criteria.  It  is  demonstrated 
that  the  internal  resistance  and  OCV  can  be  accurately  extracted 
from  the  estimated  parameters  once  the  estimated  state  error  tends 
to  zero.  Flence,  the  PE  condition  demanded  to  assure  the  asymptot¬ 
ical  stability  of  the  parameters  in  the  adaptive  control  system  can 
be  relaxed  so  that  the  convergent  performance  of  the  estimated 
internal  resistance  and  OCV  is  enhanced.  Particularly,  two  filters  are 
conceived  to  mitigate  measurement  noise  as  well  as  disturbance  so 
as  to  improve  the  reliability.  The  simulations  and  experiments  are 
both  performed  to  demonstrate  that  the  accurate  and  robust  esti¬ 
mations  of  the  internal  resistance  and  OCV  can  be  achieved.  Due  to 
the  simplicity  of  the  proposed  scheme,  it  can  be  easily  implemented 
via  electronic  circuit  design  or  embedded  system  programming. 


References 

[1  ]  B.F.W.H.  Heyer,  US  Patent  2,225,051  (1938). 

[2]  E.P.  Finger,  E.M.  Marwell,  US  Patent  4,012,681  (1975). 

[3]  T.  Kikuoka,  H.  Yamamoto,  N.  Sasaki,  K.  Wakui,  K.  Murakami,  K.  Ohnishi,  G. 
Kawamura,  H.  Noguchi,  F.  Ukigaya,  US  Patent  4,377,787  (1980). 

[4]  G.R.  Seyfang,  US  Patent  4,949,046  (1988). 

[5]  C.C.  Christianson,  R.F.  Bourke,  US  Patent  3,946,299  (1975). 

[6]  R.L.  Eby,  US  Patent  4,180,770  (1978). 

[7]  E.  Peled,  H.  Yamin,  I.  Reshed,  D.  Kelrich,  S.  Rozen,  US  Patent  4,725,784  (1984). 

[8]  K.  Muramatsu,  US  Patent  4,678,998  (1985). 

[9]  A.J.  Salkind,  C.  Fennie,  P.  Singh,  T.  Atwater,  D.E.  Reisner,  Journal  of  Power  Sources 
80(1999)293-300. 

[10]  O.  Gerard,  J.N.  Patillon,  F.  d’Alche-Buc,  Lecture  Notes  in  Computer  Science  1327 
(1997)  1095-1100. 

[11]  S.  Grewal,  D.A.  Grant,  IEEE  International  Conference  on  Telecommunications 
Energy  484(2001)  113-120. 

[12]  G.  Plett,  Journal  of  Power  Sources  134(2004)277-292. 

[  1 3  ]  J.  Lee,  O.  Nam,  B.H.  Cho,  Journal  of  Power  Sources  1 74  (2007)  9-15. 

[14]  Y.H.  Sun,  H.L.  Jou,  J.C.  Wu,  IEEE  International  Conference  on  Power  Electronics 
and  Drive  Systems  (2007)  262-266. 

[15]  P.E.  Pascope,  A.H.  Anbuky,  IEEE  Transactions  on  Power  Electronics  19(6)  (2004) 
1515-1522. 

[16]  R.J.  Grube,  Master’s  Thesis,  Department  of  Electrical  Engineering,  Wright  State 
University,  2008. 

[17]  H.  Blanke,  O.  Bohlen,  S.  Buller,  R.W.  De  Doncker,  B.  Fricke,  A.  Hammouche,  D. 
Linzen,  M.  Thele,  D.U.  Saucer,  Journal  of  Power  Sources  144  (2005)  418-425. 

[18]  B.  Hariprakash,  S.K.  Martha,  A.  Jaikumar,  A.K.  Shukla,  Journal  of  Power  Sources 
137(2004)128-133. 

[19]  B.Y.  Liaw,  R.G.  Jungst,  G.  Nagasubramanian,  H.L.  Case,  D.H.  Doughty,  Journal  of 
Power  Sources  140  (2005)  157-161. 

[20]  F.  Huet,  Journal  of  Power  Sources  70  (1998)  59-69. 

[21  ]  C.R.  Gould,  C.M.  Bingham,  D.A.  Stone,  P.  Bentley,  IEEE  International  Conference 
on  Power  Electronics  Specialists  (2008)  4381-4385. 

[22]  I.  Snihir,  W.  Rey,  E.  Verbitskiy,  A.  Belfadhel-Ayeb,  P.H.L.  Notten,  Journal  of  Power 
Sources  159  (2)  (2006)  1484-1487. 

[23]  Y.H.  Chiang,  W.Y.  Sean,  IEEE  International  Conference  on  Power  Electronics  and 
Intelligent  Transportation  System  1  (2009)  110-115. 


3932 


Y.-H.  Chiang  et  al.  /  Journal  of  Power  Sources  196  (201 1 )  3921-3932 


[24]  D.W.  Dennis,  V.S.  Battaglia,  A.  Belanger,  Journal  of  Power  Sources  1 10  (2)  (2002) 
310-320. 

[25]  J.  Newman,  K.E.  Thomas,  H.  Hafezi,  D.R.  Wheeler,  Journal  of  Power  Sources 
119-121  (2003)838-843. 

[26]  L.  Benini,  G.  Castelli,  A.  macci,  E.  Macci,  M.  Poncino,  R.  Scarsi,  IEEE  Transactions 
on  VLSI  Systems  9(5)  (2001 )  630-640. 

[27]  L.  Gao,  S.  Liu,  R.A.  Dougal,  IEEE  Transactions  on  Components  and  Packaging 
Technologies  25  (3)  (2002)  495-505. 

[28]  R.C.  Kroeze,  P.T.  Krein,  IEEE  International  Conference,  2008,  pp.  1336-1342. 


[29]  M.  Krstic,  I.  Kanellakopoulos,  P.  Kokotovic,  Nonlinear  and  Adaptive  Control 
Design,  John  Wiley  &  Sons,  1995. 

[30]  H.K.  Khalil,  Nonlinear  Systems,  Prentice  Hall,  1996. 

[31  ]  K.S.  Narendra,  A.M.  Annaswamy,  International  Journal  of  Control  45  (1 )  (1987) 
127-160. 

[32]  W.Y.  Sean,  Y.H.  Chiang,  C.H.  Wu,  Y.C.  Liang,  J.C.  Ke,  S.M.  Lo,  The  25th  Interna¬ 
tional  Battery,  Hybrid  and  Fuel  Cell  Electric  Vehicle  Symposium  &  Exhibition, 
2010. 


